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Abstract 

Results of ASCA X-ray study of central regions of medium-richness clusters of galaxies are summarized, 
emphasizing differences between cD and non-cD clusters. The intra-cluster medium (ICM) is likely to 
consist of two (hot and cool) phases within ^ 100 kpc of a cD galaxy, where the ICM metallicity is also 
enhanced. In contrast, the ICM in non-cD clusters appears to be isothermal with little metallicity gradient 
right to the center. The gravitational potential exhibits a hierarchical nesting around cD galaxies, while 
a total mass-density profile with a central cusp is indicated for a non-cD cluster Abell 1060. The iron- 
mass-to-light ratio of the ICM decreases toward the center in both types of clusters, although it is radially 
constant in peripheral regions. The silicon-to-iron abundance ratio in the ICM increases with the cluster 
richness, but remains close to the solar ratio around cD galaxies. These overall results are interpreted 
without appealing to the popular cooling-flow hypothesis. Instead, an emphasis is put on the halo-in-halo 
structure formed around cD galaxies. 
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1. INTRODUCTION 

Central regions of galaxy clusters are of particular in- 
terest with respect to the formation and evolution of 
galaxies and clusters. At or near the center of some 
clusters, we often find a single dominant elliptical galaxy 
(cD galaxy). With large extended stellar envelopes (e.g. 
Schombert 1986; Johnstone et al. 1991), these objects 
form by far the most luminous class of galaxies. Ex- 
tensive theoretical and numerical studies generally sug- 
gest that they were formed relatively early in the course 
of cluster evolution (e.g., Merritt 1985; Dubinski 1998; 
Ghigna et al. 2000), but detailed comparison with the 
observation is yet to be performed. In other clusters, 
we find several giant elliptical galaxies instead of a cD 
galaxy. The origin of the difference between clusters 
with and without cD galaxies is not yet clear. 



As noticed by Jones, Forman (1984), the presence of 
a cD galaxy in a cluster also affects the properties of 
the intra-cluster medium (ICM), which fills the intra- 
cluster space and emit intense X-rays (Sarazin 1988). 
Around a cD galaxy, the soft X-ray surface brightness 
often exhibits a strong excess above a (3 model fitted to 
the outer-region brightness profile (e.g. Jones, Forman 
1984; Edge et al. 1992); we may call the phenomenon 
"central excess emission" (CEE). In addition, the X-ray 
emission from such an environment often involves a spec- 
troscopic cool component, hereafter called central cool 
component (CCC), emitted from plasmas with a repre- 
sentative temperature of T 1 keV (e.g. Fabian 1994; 
Fabian et al. 1994; Matsumoto et al. 1996). So far, the 
CEE that is a spatial effect and the CCC that is a spec- 
troscopic one have been interpreted as different aspects 
of the same phenomenon, i.e., radiative cooling process 
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and associated continuous ICM in-flow, termed cooling 
flow (CF; e.g. Fabian 1994), proceeding in the densest 
part of the ICM. 

The CF scenario was motivated by the short radiative 
cooling times of the ICM inferred from X-ray observa- 
tions, and reinforced by the CCC and CEE phenom- 
ena. Theoretical treatment of CFs has been developed 
to such a degree of sophistication as inhomogeneous CFs 
(Nulsen 1986; Fabian 1994) and isothermal CFs (Nulsen 
1998). However, the CF hypothesis is still subject to sev- 
eral basic problems. The heat conduction must somehow 
be suppressed in order for a CF to develop (Takahara, 
Takahara 1979). The fate of the large amount of cooled 
gas is controversial: the possible X-ray absorbing gas 
(e.g. White et al 1991) may not be massive enough to 
act as a full sink for CFs (Allen, Fabian 1997), and the 
presence of such an absorber itself is sometimes ques- 
tioned (Sarazin 1997: Huang, Sarazin 1998; McNamara 
et al. 2000). The searches for neutral hydrogen (e.g., 
Dwarakanath et al. 1994) and coronal emission lines 
(e.g., Yan, Cohen 1995) have been unsuccessful. Fur- 
thermore, numerical calculations of the cluster evolu- 
tion predict too high X-ray luminosities when the ICM 
cooling is fully taken into account (Suginohara, Ostriker 
1998). Although Cen, Ostriker (1999) predict that a 
major fraction of baryons in the universe is in the form 
of warm gas of temperature 10^~^ K, this component 
is predicted to diffusely distribute in the inter-galactic 
space rather than forming a cooling condensation. 

X-ray observations by ASCA (Tanaka et al. 1994), 
with its imaging capability in the previously unexplored 
energy range up to 10 keV and its high spectroscopic 
performance, have renewed our understanding of the 
CCC phenomenon. Although we can clearly resolve the 
CCC in the ASCA spectra taken from central regions of 
many cD clusters (e.g. Fabian et al. 1994), the strongest 
CCC has been found to account for only a minor frac- 
tion of the X-ray emission therein (subsections 2.1 and 
2.2). This implies that the effect of ICM cooling was 
considerably overestimated previously (subsection 2.5). 
Furthermore, the ICM around cD galaxies often exhibits 
a marked metallicity increase and characteristic chemi- 
cal compositions (subsection 4.3), neither of which fits 
into the CF hypothesis. Accordingly, we can more nat- 
urally interpret the CCC as an X-ray emitting inter- 
stellar medium (ISM) associated with the cD galaxy 
(subsection 2.6), rather than as the cooling portion of 
the ICM. This in turn allows us to construct a realis- 
tic scenario of production, confinement, and transport 
of heavy elements, from the scales of individual galaxies 
to the cluster-wide scale. 

Spatial studies of nearby clusters with ASCA have also 
yielded a surprise; the CEE that was known previously 
in soft X-rays has also been detected in hard X-rays up to 
^ 10 keV (subsection 3.1). The revealed CEE properties 
in fact depend little on the X-ray energy. We accordingly 
regard the CEE as a manifestation of particular shapes 
of the gravitational potential at the cluster center (sub- 
sections 3.3 and 3.4), rather than the ongoing ICM cool- 



ing. This allows us for the first time to compare in close 
detail the X-ray determined potential shapes with those 
suggested by A''-body simulations (Navarro et al. 1996; 
Fukushige, Makino 1997; Moore et al. 1998; Ghigna et 
al. 2000), or those derived from gravitational lensing 
measurements (e.g. Wu et al. 1998). 

The present paper is meant to provide an overall 
summary of these new ASCA results, mainly based on 
six PhD thesis works (Ikebe 1995; Matsushita 1997; 
Fukazawa 1997; Xu 1998; Tamura 1998; Ezawa 1998). 
Although many of these results have already been pub- 
lished individually, here we assemble them together as 
building blocks, and attempt to construct a novel com- 
prehensive view of the physics in the cluster core regions. 
We devote Sections 2, 3, and 4 to the descriptions of 
the ICM temperature structure, the spatial properties of 
the ICM, and the ICM metallicity, respectively. There, 
we incorporate brief comparison with the early Chandra 
and XMM-Newton results. These results are combined 
together in Section 5. 

In order to avoid substructures often seen in rich 
clusters, we mainly investigate relatively poor clusters, 
with ICM temperature below ~ 6 keV. Unless otherwise 
stated, we use the 90% confidence error regions through- 
out, and assume the cluster to be spherically symmetric. 
For the sake of easy comparison with previous results, 
we assume the Hubble constant to be Hq = 50 h^o km 
Mpc~^, even though it somewhat differs from the 
most recent determinations. We also employ the solar 
abundance ratios from Anders, Grevesse (1989), with 
Fe/H= 4.68 x 10~^ by number. We denote the three- 
dimensional radius R, while the projected radius r. 

2. PROPERTIES OF THE CENTRAL COOL 
COMPONENT 

The CCC properties were so far studied in a lim- 
ited energy band below ~ 3 keV. In addition, few such 
observations were performed with high spectroscopic 
resolution, except those with the Einstein FPCS (e.g., 
Canizares et al. 1982), the Einstein SSS (e.g., Rothen- 
flug et al. 1984), and the BBXRT (e.g., MacKenzie et 
al. 1996). These limitations, coupled with projection 
effects, made it difficult to unambiguously quantify how 
the ICM temperature decrease actually takes place (e.g., 
Thomas et al. 1987). The superior energy resolution of 
ASCA, together with its imaging capability up to ~ 10 
keV, have therefore enabled us to acquire novel informa- 
tion on the CCC phenomenon. 

2.1. The Centauru.s cluster 

The Centaurus cluster, with a relatively good circular 
symmetry, a low redshift (0.011), and the most promi- 
nent CCC among nearby clusters (Matilsky et al. 1985; 
Allen, Fabian 1994; Fabian et al. 1994; Fukazawa et 
al. 1994; Ikebe 1995; Ikebe et al. 1999), provides an 
ideal opportunity to investigate the ICM temperature 
structure. According to the ROSAT data (Allen, Fabian 
1994), the CCC of this cluster is seen within a projected 
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radius of r ~ 5' (~ 100 h^Q kpc). Fabian et al. (1994) 
described the ASCA results on the CCC in terms of the 

CF scenario. 

The same ASCA data, taken with the GIS (Gas Imag- 
ing Spectrometer; Ohashi et al. 1996; Makishima et 
al. 1996) and the SIS (Solid-State Imaging Spectrom- 
eter; Burke et al. 1994; Yamashita et al. 1999), were 
also analyzed by Fukazawa et al. (1994), Ikebe (1995), 
and Ikcbc ct al. (1999), with increasing degree of so- 
phistication. Ikebe et al. (1999) also incorporated the 
ROSAT PSPC data. These authors have established 
"two-temperature (2T)" picture of the ICM of this fasci- 
nating cluster: the spectra accumulated over concentric 
annular regions around the cluster center can be fitted 
by a sum of cool and hot emission components, with a 
temperature of Tc ^ 1 keV and Th ^ 4 keV, respectively. 
Since the relative contribution of the cool component 
clearly diminishes for r > 7', tlic^ cool component can be 
identified with the CCC. In addition, these authors have 
obtained the following two important new findings. 

One is that the spectra taken from the central (r < 5') 
region contains not only the cool component, but also a 
significant amount of hot component, which may have 
escaped detections by the low-energy imaging instru- 
ments. Fukazawa et al. (1994) argued that the ob- 
served hot emission indeed comes mostly from the three- 
dimensional cluster core region rather than from the 
foreground or background off-center regions along the 
line of sight. Therefore, the apparent 2T property does 
not result from projections of a single- phase ICM (i.e., 
the temperature being a single-valued function of R) 
having a central temperature decrease, but from multi- 
phase nature of the ICM in the core region of this clus- 
ter. This view was reinforced by Ikebe (1995), who 
showed that single-phase models with central temper- 
ature drops generally undcrpredict the hard X-ray flux 
from the cluster center region. Convincing support has 
been obtained by Ikebe et al. (1999); the ICM within 
a three-dimensional radius i? ~ 1.'5 (^ 30 h'^Q kpc) 
cannot be isothermal, and at least two components of 
different temperatures must be involved there. 

The other result concerns the number of different tem- 
perature components in the cluster core region. The two 
temperatures from the 2T fit, and Th, are both vir- 
tually constant as a function of r (Fukazawa ct al. 1994; 
Makishima 1994a), which would not generally occur if 
many components with different temperatures and dif- 
ferent angular distributions were involved. Furthermore, 
the ASCA and ROSAT spectra from the cluster center 
have been fitted well {x^/v = 119/121) with the 2T 
model (Ikebe ct al. 1999), while the fit became imac- 
ceptable jv = 159/122) when the cool component 
was replaced by the CF spectral model. The conclusion 
remained unchanged even by applying a separate free ex- 
cess absorption to the CF component, which gave only a 
slightly improved fit {y^ /y = 142/121) together with a 
best-fit excess column of ^ 1.3 x 10^^ cm^^. Therefore, 
there seems to be only two major temperature compo- 
nents, rather than a continuous distribution of emission 



measure as a function of temperature. 

We hence regard the 2T modeling of Centaurus as 
physically meaningful, rather than just a convenient de- 
scription of the data. Actually, the best-fit 2T solution 
obtained by Ikebe et al. (1999) can consistently and 
simultaneously describe the spatial and spectroscopic 
properties of the existing ASCA GIS/SIS and ROSAT 
PSPC data of this cluster. Specifically, it dictates the 
following features. 

1. The two temperatures are Th = 3.9 ± 0.1 keV and 
Tc = 1.4 ±0.2. 

2. The cool phase, or CCC, is localized to within R ~ 
5' (~ 100 H'^q ) of the cluster center, and exhibits a 
0.5-4 keV luminosity of 1.0 x 10^^ h^^ ergs s'^. 

3. Assuming a pressure balance between the two 
phases, the CCC is calculated to have a mass of 
(4 — 5) X 10^*^ Mq, and an average volume filling 
factor of 0.06-0.08 within 60 h^^ kpc. Thus, the 
central region is mainly occupied by the hot phase. 

4. Toward the center, the ICM metallicity increases 
significantly in both phases, up to 1 — 1.5 solar abun- 
dances. The central excess iron contained in the two 
phases amounts to ^ 1 x 10^ Af©. 

5. There is no evidence of excess X-ray absorption, 
beyond typical upper limits of ~ 1 x 10^^ cm~^. 

6. The gravitational potential becomes deeper at the 
center than a King-type potential, showing either a 
hierarchical structure or a central cusp. 

The first four items reconfirm or update the results ob- 
tained previously by Fukazawa et al. (1994) and Ikebe 
(1995), while the second item agrees with Allen, Fabian 
(1994). 

Thus, the first significant result of the present paper 
is the successful review of the detailed 2T model for the 
Centaurus cluster, with an implication that the clus- 
ter core volume is mostly filled with the hot phase even 
though there exists a significant CCC. 

The spatial co-existence of two distinct plasma com- 
ponents is also found in the case of diffuse Galactic X-ray 
emission, called Galactic ridge X-ray emission (GRXE; 
Koyama et al. 1986; Kaneda et al. 1997). Using 
the ASCA Galactic-plane survey data, Kaneda (1997) 
analyzed the GRXE surface brightness variation along 
the Galactic longitude. Employing cores-correlation and 
fractal dimension analyses, he has unambiguously shown 
that the 0.8-10 keV GRXE comprises two emission com- 
ponents, with characteristic temperatures of ~ 0.9 keV 
and 3-7 keV. 

2.2. Two-temperature properties of a, sample of clusters 
In clusters other than Centaurus, the CCC is generally 
not strong enough to unambiguously discriminate the 2T 
picture from other interpretations, e.g., single-phased 
ICM with a temperature gradient, or a highly multi- 
phased ICM with the temperature distributing over a 
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wide range. Nevertheless, in the case of the Centaurus 
cluster (subsection 2.1), the 2T model has been shown 
to be more appropriate than these alternatives. Since 
there is no particular reason to regard Centaurus as ex- 
ceptional, we again adopt the 2T modeling as a working 
tool, and examine whether it can describe the ASCA 
data of other clusters. 

Actually, the ASCA SIS spectrum from a central re- 
gion of the Perseus cluster, with a strong CCC, was re- 
produced reasonably well with the 2T model, although 
the "hot-phase plus CF" model was similarly successful 
(Fabian et al. 1994). The 2T formalism has also been 
successful for central regions of the Hydra-A cluster and 
Abell 1795, as reported by Ikebe et al. (1997b) and Xu 
et al. (1998) respectively. The X-ray emission from M87, 
at the center of the Virgo cluster, has also been mod- 
eled successfully by Matsumoto et al. (1996) in the same 
way, using Th = 2.4 - 3.2 keV and = 1.1 - 1.5 keV, 
although the narrower separation between the two tem- 
peratures makes it rather difficult in this case to assess 
the uniqueness of the 2T formalism. We here remem- 
ber the conclusion derived previously by Canizares et 
al. (1982) through an elaborate analysis of the Einstein 
FPCS data, that the X-ray emitting plasma in M87 is 
highly multi-phased, with the temperature distributed 
from ^ 0.3 keV to ~ 3 keV or higher. However, their 
conclusion is heavily dependent on the Fe-L complex, of 
which the atomic model calculations later turned out to 
be problematic (e.g., Fabian et al. 1994; Matsushita et 
al. 1997). Accordingly, reservations should be put on 
the interpretation by Canizares et al. (1982). 

In addition to these individual attempts, Fukazawa 
(1997) and Fukazawa et al. (1998, 2000) took a statis- 
tical approach. These authors used a sample of 40 clus- 
ters observed with ASCA, and accumulated X-ray spec- 
tra for each object separately over the central (within 
r 100 h'^Q kpc) and outer regions. The sample defi- 
nition and details of the data analysis are described in 
Fukazawa et al. (1998, 2000). For all clusters, the outer- 
region spectra were fitted successfully with a single- 
temperature plasma emission model. In contrast, the 
same modeling failed for the central-region spectra from 
about half of the objects in the sample. These spectral 
fits have been improved significantly by employing the 
2T formalism, in which the value of Th is set equal to 
the outer-region ICM temperature of the same cluster. 
These clusters hence have statistically significant CCC, 
for which the 2T formalism gives a successful descrip- 
tion. The 2T parameters derived by Fukazawa (1997) 
and Fukazawa et al. (2000) are quoted in table 1, after 
discarding too rich (Th > 6.0 keV) or too poor (Th < 2.0 
keV) clusters: the former criterion is to avoid objects 
with substructures. Our reduced sample comprises 20 
objects, including, e.g., Centaurus, Virgo, Hydra-A, and 
Abell 1795, but not Perseus which has Th > 6 keV. Thus, 
the CCC temperature is found at Tc = 1.1 — 2.2 keV, in 
agreement with the results on the Centaurus and Virgo 
clusters. These results extend the validity of the 2T pic- 
ture to a larger sample. 



Also listed in table 1 is the ratio Qc/Qh, where Qc 
and Qh sue emission integrals of the cool and hot com- 
ponents, respectively, both integrated over the central 
^ 100 h^Q kpc. These are related to the hot-phase and 
cool-phase luminosities, Th(Ai?) and L\^{AE) respec- 
tively, emergent from that region in a specified energy 
band AE, as 

Lh(AS) = Qh-Ah(AS) , L,{AE) = gc-Ae(A^) , (1) 

where Ah(A£:) = A(Th,Zh;A£') and Ac(A£') = 
A(Tc, Zc, AE) are the band-limited cooling functions for 
the hot and cool phases, respectively, with and Zc 
denoting their metallicities. 

Assuming the CCC to be confined within a volume V 
with an average volume filling factor 77 (Fukazawa et al. 
1994; Makishima 1994b), the emission integrals can be 
written as 



Qh = (1 - 7?)V , Qc^{nl)fiV 



(2) 



where Uh and Uc are plasma densities of the hot and cool 
phases respectively, and the bracket means the spatial 
average over V. Employing the assumption of pressure 
balance between the two phases (Fukazawa et al. 1994; 
Makishima 1994b) as 



ricTc = 7ihTh , 
equation (0) yields 



V = 



1 



Qc 



(3) 



(4) 



[the same as eq.(9) of Ikebe et al. (1999)]. 

Taking the Centaurus cluster for example, the values 
listed in table 1, Qc/Qh ~ 1-0 within r = 100 h^Q kpc 
and Th/Tc ~ 2.7, give 77 ~ 0.1 via equation (|). This 
roughly agrees with the results from the more detailed 
analysis by Ikebe et al. (1999) quoted in subsection 2.1. 
Similarly, in table 1, we give the value of 77 for each 
cluster calculated via equation (||). Thus, rj takes rather 
small values, falling below 0.1 except for Abell 496, Cen- 
taurus, 3A 0335-1-096, and Virgo. This implies that the 
core regions of clusters in table 1 are occupied predom- 
inantly by the hot phase even if there exists a CCC. 
Hereafter, we employ the 2T formalism as our standard 
tool, and identify the CCC with the cool component. 

From equations (|l|), (||), and (^, we immediately ob- 
tain 



L,{AE) _(Q. 



Lh(AS) \Q 



Ae 

Ah 



Ae 

Ah 



1 



.(5) 



This allows us to relate Tc/Th, Qc/Qh, and rj to one 
another. Again taking Centaurus for example, substi- 
tution of 7y 0.1 into equation (^ yields Lc/L^ ~ 
0.8Ac/Ah, and actual calculation of the cooling function 
(using Th = 3.8 keV and Tc = 1.3 kcV, and assuming 
Zh ~ 1) gives the Tc/^h ratio of 1.1, 0.16, and 

0.68, in the 0.5-3, 3-10, and 0.5-10 keV bands, respec- 
tively; these values are consistent with those calculated 
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by Ikebe et al. (1999) for the Centaurus cluster. Thus, 
in soft energies where Ac > Ah, the CCC luminosity Lc 
of a cD cluster can be comparable to the hot-phase lu- 
minosity Lh from the cluster core region. In contrast, 
in higher energies or in a sufficiently wide energy range, 
Lc obviously falls much below L^, because 77 << 1 and 
Ac « Ah. 

The Chandra and XMM-Newton observations are pro- 
ducing a series of interesting results on the ICM temper- 
ature structure, that there is little X-ray emission com- 
ponent with the temperature below a certain lower cutoff 
value that differs from object to object. The cutoff tem- 
perature is reported to be 2.7 keV for Abell 1835 (Peter- 
son et al. 2001), 1.5 keV for Sersic 159-03 (Kaastra et al. 
2001), and 2.5 keV for AbeU 1795 (Fabian et al. 2001; 
Tamura et al. 2001) which is close to our Tc (2.04 keV; 
table 1). These results suggest that the X-ray spectra 
of each cluster are characterized by two temperatures, 
i.e., the lower cutoff value and the outer-region averaged 
ICM temperature. Thus, our 2T interpretation appears 
to be valid for these new data as well, although it is yet 
to be examined whether or not these new missions are 
detecting emission components with the temperature in 
between the two characteristic values. 

2.3. The central cool component and cD galaxies 

Although found in many clusters, the CCC is subject 
to a large scatter, from object to object, in its promi- 
nence relative to the whole cluster emission. For ex- 
ample, compared to the Centaurus cluster, the nearby 
cluster Abell 1060 is rather similar in distance, optical 
richness, ICM temperature, and the overall X-ray lu- 
minosity. Nevertheless, its ICM is quite isothermal at 
3.1 ±0.2 keV right to the center, with little evidence for 
CCC. More quantitatively, Tamura et al. (1996, 2000) 
and Tamura (1998) constrained the 0.3-5 keV luminos- 
ity of CCC in this cluster to be < 6 x 10^^ h^^ ergs s~^ 
over a central region of r < 3' (< 60 h^^ kpc). This 
falls by more than an order of magnitude below that 
of the Centaurus cluster. The same conclusion can be 
derived by comparing the two clusters in terms of the 
Qc/Qh ratio given in table 1; it is < 0.1 for Abell 1060, 
while ~ 1 for Centaurus. In addition, Abell 1060 has 
long been know as a typical weak-CEE object (Jones, 
Forman 1984; see subsection 3.4). 

What makes the two clusters so different? Although 
CCC could be disrupted by mergers (e.g. Fabian 1994), 
Abell 1060 is even more relaxed and rounded than Cen- 
taurus, with least evidence for recent mergers. Several 
authors, including Tamura et al. (1996), instead ascribe 
the difference to the properties of their central galax- 
ies. The Centaurus cluster is classified as Bautz-Morgan 
(B-M; Bautz, Morgan 1970; Abell et al. 1989) type I- 
II, and hosts a single giant galaxy, NGC 4696, which 
completely dominates the central cluster region up to 
~ 200 kpc in radius. In contrast, Abell 1060 has 
B-M type III, and host two central elliptical galaxies, 
NGC 3309 and NGC 3311, with a projected separation 
of only ~ 35 h'^Q kpc. There is yet another bright spiral 



galaxy, NGC 3314, at a separation of ~ 100 h^Q kpc 
from NGC 3311. None of the three dominates the cen- 
tral region of Abell 1060. Tamura et al. (1996) argue 
that these morphological properties are responsible for 
the X-ray difference between the two clusters; the idea 
may be traced back to Jones, Forman (1984). 

In order to investigate the above suggestion, we here 
subdivide the objects in table 1 into two subsamples, cD 
and non-cD clusters. In order to avoid tautology, we uti- 
lize solely the optical information: we define cD clusters 
as those of B-M type I, I-II, or II, while non-cD ones as 
those with B-M type II-III or III. The only exception is 
the Virgo cluster; it is of B-M type III, but we classify 
it as a cD cluster, because M87 is widely regarded as 
its cD galaxy. Although the B-M type is unavailable for 
3A 0335-1-096 and the Hydra-A cluster, we classify them 
as cD clusters based on their Rood-Sastry classification 
as "cD" (Rood, Sastry 1971; Struble, Rood 1987). Our 
classification is not much different from that of Fukazawa 
(1997) and Fukazawa et al. (2000), who mainly utilized 
the Rood-Sastry morphology. 

In figure |l], we plot the Qc/Qh ratio taken from ta- 
ble 1 as a function of Th. It is thus clear that the cD 
clusters exhibit systematically higher Qc/Qh ratios than 
the non-cD ones. We hence conclude that the CCC is 
selectively seen among cD clusters. Although the gen- 
eral association of CCC with cD galaxies has long been 
known (e.g., Jones, Forman 1984; Edge et al. 1992; 
Fabian 1994), figure |l] allows us for the first time to 
quantitatively discriminate cD and non-cD clusters in 
terms of the 2T picture. 

2.4. Mass and luminosity of the central cool component 
The 2T modeling allows us to calculate the mass as- 
sociated with the CCC (i.e., the cool-phase mass), Mc, 
as performed by Fukazawa et al. (1994) for Centaurus. 
However for this purpose, we need to quantify the ICM 
density profile. For this reason, we cannot derive Mc 
for all objects in table 1. Instead, we refer to individ- 
ual publications dealing with selected cD clusters, and 
compile the reported values of Mc in figure |2] against 
Lc. Thus, the values of Mc of the selected cD clusters 
are at most several times 10^'^ Mq, or only a few to ten 
percent of the stellar mass in a typical cD galaxy. 

In figure ^ we also show masses and X-ray lumi- 
nosities of the ISM (inter-stellar medium; Forman et al. 
1985; Trinchieri, Fabbiano 1985; Canizares et al. 1987) 
associated with non-cD elliptical galaxies, obtained by 
Matsushita (1997) using ASCA. Thus, the CCC masses 
are only a few times higher than the ISM masses of the 
non-cD elliptical galaxies with the highest X-ray lumi- 
nosities (e.g. NGC 4472, NGC 4636, NGC 5846, and 
IC 4296). The CCC and the ISM emission are similar to 
each other in temperature and angular extent as well. 
Furthermore, the CCC is predominantly seen around 
cD galaxies (subsection 2.2). These comparisons reveal 
a close similarity between the CCC in cD clusters and 
the ISM emission from X-ray luminous elliptical galaxies 
(see subsection 2.6 for a further discussion). Neverthe- 
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less, in figure ^ the CCC exhibits significantly higher 
luminosities than the ISM of ellipticals. 

From equations (Q) and (||), the cool-phase luminosity 
scales as 



(6) 



where definitions of the symbols are the same as in equa- 
tions @) and (|). Therefore, even M^, V, T^, and Z (the 
latter two affecting Ac) are kept constant, the luminos- 
ity can be increased by decreasing ij. The behavior of 
CCC in figure |2] can be explained in this way, because 
the 2T interpretation implies that the CCC actually has 
very small values of 77 (subsection 2.2). Presumably, the 
CCC forms blobs or filaments around each cD galaxy, 
compressed by the ambient hot ICM. 

We thus conclude that the CCC of a cD galaxy has a 
close resemblance to the ISM of X-ray luminous ellipti- 
cal galaxies with respect to the mass and spatial extent, 
but has a higher luminosity presumably because it is 
confined to achieve a lower filling factor. The small val- 
ues of Mc and the relatively high values Lc lead to rather 
short radiative cooling times for the cool component, as 
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A similar story applies to another well known cooling- 
flow cluster, Abell 1795. With Th = 6 keV, this cluster is 
reported to have M = 478 Mq yr~i with EXOSAT 
(Edge et al. 1992), or M - 500 h^^ Mq yr"! with 
ROSAT (Briel, Henry 1996; Allen, Fabian 1997). In con- 
trast, the ASCA data yielded = (1.4±0.4) x lO''* h^^ 
erg s-i in the 0.5-3 keV band (Xu et al. 1998), which 



Mq yr via bolometric cor- 
Consistently, the CF-model 



2.5. Amount of cool emission 

The cool-component luminosities measured with 
ASCA, summarized in table 1 and figure ^, are con- 
siderably lower on average than was thought previously. 
To compare our results with previous ones, let us convert 
the bolometric cool-component luminosity Lj?°' mea- 
sured with ASCA to the mass deposition rate M in terms 
of the CF scenario, utilizing a theoretical relation of 



(Fabian 1994) where k is the Boltzmann constant, Wp 
is the proton mass, and /x ~ 0.6 is the mean molecular 
weight. The ICM is assumed to start cooling from Th- 

The Hydra-A cluster is one of the prototypical cooling- 
flow clusters with the reported mass deposition rate 
reaching M = 
1992), - 600 



Mq yr-^ (Edge et al. 
Mq yr-i (David et al. 1990), 
or - 270 Mq yi'^ (Allen, Fabian 1997; using 

ROSAT). However, a joint 2T fit to the ASCA and 
ROSAT spectra of the Hydra-A cluster yielded a rather 
weak cool component with Lc ~ 5 x 10^'^ h^Q erg cm~^ 
s-i in the 0.5-3 keV band (Ikebe et al. 1997b). Then, 
from equation (^) and ^ A keV, and after bolo- 
metric correction, we obtain M ~ 60 h'^Q Mq yT~^, 
which falls by 5-10 times below the previous estimates. 
To examine if this result depends on our 2T model- 
ing, Ikebe et al. (1997b) further fitted the same spec- 
tra directly with the CF spectral model of Mushotzky, 
Szymkowiak (1988), and obtained a similarly low value 
of M = (60 ± 30) h~Q Moyr-i. Applying a separate 
excess absorption to the CF component did not change 
the result significantly. 



translates to M ^ 150 
rection and equation (H) 
fit to the ASCA SIS spectra yielded M - 131 h^Q Mq 
yr~^ (Fabian et al. 1994). Thus, the values of M de- 
rived with ASCA are again 3^4 times lower than the 
previous estimates. 

Yet a third interesting example is Abell 1060 discussed 
in subsection 2.3. This cluster was thought to host a 
rather weak cooling flow with Al = (2.4 — 15) h^Q Mq 
yr-i (Singh et al. 1988), or M = 6 h^^ Mq yr^^ (Stew- 
art et al. 1984). However, the ASCA upper limit on the 
CCC (tabic 1) is so low that it gives, after bolometric 
correction, M < 1 h^^ yr~^ in terms of equation (||). 
Again, this falls by an order of magnitude below the 
previous estimates. 

Including these particular examples, we compare in 
^y^gurc ^ the mass deposition rates determined spec- 
troscopically with ASCA, against those derived in soft 
X-rays mainly employing the surface brightness pro- 
files. Although the two estimates agree on the clus- 
ters with the strongest CCC (Centaurus, Virgo, and 
3 A 0335-1-096), for other clusters the ASCA values fall 
systematically and significantly below the previous ones. 
This indicates that the past imaging spectroscopy in 
soft X-rays (especially, the commonly used deprojec- 
tion analysis) has considerably overestimated the cooling 
mass deposition rates of clusters. 

This conclusion is being reconfirmed by a series of 
new results from Chandra and XMM-Newton. For ex- 
ample, analysis of the Chandra data of the Hydra- 
A cluster has led McNamara et al. (2000) to derive 
M = (34 ± 5) Mqji-^ within i? = 74 kpc, for Hq = 70 
km s~^ Mpc~^: the value would not greatly exceed 
100 M0yr~^ even if we integrate up to ~ 200 kpc. 
The spectra do not require excess absorption above the 
Galactic column of 2 x 10^° cm~^. These results have 
further been detailed by David et al. (2001). Similar 
deficits of the cool plasma are reported on Abell 1795 
(Tamura et al. 2001), AbeU 1835 (Peterson et al. 2001), 
and Sersic 159-03 (Kaastra et al. 2001), all being typical 
CF clusters. The reported scarcity of low-temperature 
component is independent of the validity of our 2T pic- 
ture. 



2. 6. The nature of the central cool component 

In the high-quality ASCA spectra taken from central 
regions of cD clusters, we have thus successfully resolved 
the CCC as a spectroscopic cool component. The CCC, 
or the cool phase, is locaHzed around the cD galaxy (sub- 
section 2.3), presumably forming blobs or filaments im- 
mersed in the vast sea of the hot phase (subsection 2.4). 
This picture apparently agrees with the scenario of in- 
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homogeneous CF (Nulsen 1986; Fabian 1994). However, 
invoking the CF hypothesis in its original form may no 
longer be appropriate, because we have at the same time 
discovered that the ICM deposition rate due to radiative 
cooling was significantly overestimated in the past, by up 
to an order of magnitude (subsection 2.5). Then, what 
is the nature of CCC? 

Considering that non-cD elliptical galaxies in clusters 
(e.g., NGC 4406 and NGC 4472 in the Virgo cluster; 
Forman et al. 1985; Awaki et al. 1994) usually possess 
their own ISM, a cD galaxy, which is far less subject to 
ram-pressure stripping, must have its own ISM as well. 
Then, by extending the argument of subsection 2.4, it is 
most natural to interpret the central cool phase as the 
ISM associated with the cD galaxy (Makishima 1996, 
1997b, 1999b; Ikebe et al. 1999), rather than as a cooling 
portion of the ICM. Specifically, we have the following 
pieces of evidence in support of this view; (1) the CCC 
is statistically associated with cD galaxies (subsection 
2.3, figure [l]); (2) the CCC is positionally centered on 
cD galaxies (e.g., Lazzati, Chincarini 1998) even if they 
are sometimes offset from the dynamical cluster centers; 
(3) both the CCC and the ISM emission have a similar 
angular extent of several tens kpc; (4) the CCC temper- 
ature, Tc = 1 ~ 2 keV, is close to that of an elliptical's 
ISM, which in turn is consistent with the potential depth 
of giant elliptical galaxies; (5) the estimated mass of the 
cool phase is comparable to the ISM mass of ordinary 
elliptical galaxies (subsection 2.4, figure The only 
difference is the considerably higher X-ray luminosity 
of CCC compared to the ordinary ISM emission, but 
this can be explained as a result of compression by the 
hot phase (subsection 2.4). The important point is that 
the CCC has a lower temperature because it reflects the 
shallower potential depth of a cD galaxy, rather than 
because it is radiatively cooling. Later, we reinforce this 
interpretation from metallicity arguments. 

The above new picture, however, involves two imme- 
diate problems. One is theoretical; what sort of energy 
input (heating) prevents the cool phase from radiative 
cooling collapse, and how are the two phases kept ther- 
mally insulated from each other? The heating mecha- 
nism to be invoked, though not required to have as large 
a luminosity as thought previously, must be rather effi- 
cient to balance the rapid radiative cooling of equation 
. Furthermore, the heating must balance the cooling 
in a stable manner, in order to sustain the cool phase. 
These are indeed the fundamental issues which under- 
lie the CF hypothesis (Fabian 1994); we briefly discuss 
them in Section 5. 

The other problem is observational; what caused the 
previous overestimates of M? An obvious possibility is 
ignorance of the dominance of the hot phase in the clus- 
ter core regions, which in turn was due to the limited 
energy range and insufficient spectral resolution of the 
previous imaging X-ray observations. Observers tend to 
attribute all the X-ray flux measured from the core re- 
gion, which is in fact dominated by the hot component, 
to the cooling flux. However, this idea alone is inade- 



quate to fully answer the question. Actually, even using 
the ROSAT PSPC data alone, discrepant values of M 
are sometimes assigned to a single object depending on 
the data analysis method. For example, M < 80 Mq 
yr^^ is obtained through analysis of the PSPC spectra 
of Abell 4059, whereas the surface brightness analysis of 
the same PSPC data yields M ~ 184 Mq yr'^ for the 
same object (Huang, Sarazin 1998). Similarly, values 
of M as high as ~ lO'^ Mq yr~^ are reported for some 
distant clusters via deprojection analysis of the ROSAT 
data, even though the spectral evidence for CCC is poor 
(e.g., Schindler et al. 1997). Therefore, the remaining 
clue may reside in the surface brightness proflles, espe- 
cially the CEE phenomenon. This urges us to study, in 
the next Section, spatial properties of the cluster X-ray 
emission using ASCA. 

3. GRAVITATIONAL POTENTIAL SHAPE 
IN CLUSTER CENTERS 

Imaging soft X-ray observations have established the 
CEE (central excess emission) as a ubiquitous phe- 
nomenon among cD clusters. For example, analysis of 
the ROSAT X-ray surface brightness proflles for a large 
sample of clusters, using single-/? modeling, reveals a 
tight correlation between (3 and the core radii (Pownall, 
Stewart 1996; Neumann, Arnaud 1999). This indicates 
the prevalence of CEE, because both these quantities 
are known to take rather small values when a single-/? 
model is forced to fit a brightness profile with a signif- 
icant CEE (Makishima 1995). However, before ASCA, 
there has been essentially no CEE investigation at ener- 
gies above 3 keV. 

3.1. Central excess emission in hard X-rays 

Generally, the X-ray volume emissivity e of the ICM 
in a given energy band Ai? is expressed as 

e(i?; l^E) = n{Rf ■ A{T, Z; AE) . (9) 

For reference, volume integration of this equation for 
each ICM phase, taking into account the filling factor, 
together with equation (||), yields equation (|l|). Thus, e 
is directly proportional to n^, and hence the CEE phe- 
nomenon implies an excess ICM density in the cluster 
core region above the prediction of a /3 model. This 
statement remains valid even if the ICM is deviated sig- 
nificantly from isothermality, since the cooling function 
A(T, Z; AE) depends only weakly on T in the relevant 
parameter regime. 

The implied excess ICM density in the cluster core 
region, in turn, may be produced by either of the fol- 
lowing two mechanisms. One is that the ICM pressure 
distribution has a flat core without central excess, and 
the increase in n(R) is compensated by a decrease in the 
ICM temperature, such as is invoked in the CF hypoth- 
esis. The other possibility is that the ICM is relatively 
isothermal, and hence the ICM pressure itself exhibits 
an excess in the central region. So far, the attention 
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has been paid predominantly to the former mechanism, 
based on a belief that the CEE phenomenon is a spatial 
counterpart to the CCC phenomenon. Observationally, 
in energies below ~ 3 keV, these alternatives produce 
very similar effects and are difficult to distinguish. How- 
ever, in sufficiently high energies, their effects become 
distinct: the former predicts a central deficit in the hard 
X-ray brightness because the hot phase is displaced by 
the increasing cool-phase contribution, while the latter 
obviously predicts a central excess in hard X-rays. 

In order to examine the issue, Ikebe et al. (1997b) and 
Xu et al. (1998) respectively analyzed the ASCA data of 
the Hydra-A cluster and Abell 1795, both known to host 
strong soft-band CEE (e.g. Jones, Forman 1984; David 
et al. 1990). By fully taking into account the instrumen- 
tal responses of ASCA, they have discovered that these 
clusters surprisingly exhibit, even in high energies above 
~ 3 keV, a clear CEE above the (3 models describing 
their outer-region X-ray profiles. Furthermore, in both 
objects, the relative prominence of CEE has been found 
to be roughly energy independent, from the ROSAT 
range (0.2-2 keV) up to - 10 kcV (figure 9 of Ikebe 
et al. 1997b and table 3 of Xu et al. 1998). According 
to Xu (1998), the same is approximately true of AWM7 
and Abell 2199, which are also known to show promi- 
nent CEE in soft X-rays. Since the CCC contribution in 
these objects is negligible above ~ 3 keV, the approxi- 
mately color-independent CEE unambiguously indicates 
the presence of central excess pressure, supporting the 
idea that the ICM pressure exhibits an increase at small 
radii rather than that the ISM pressure has a flat core 
(and hence the gas temperature decreases). 

Then, what happens in clusters with very strong CCC, 
such as Centaurus, in which the soft-band CEE is also 
strong (Matilsky et al. 1985)? Using ASCA, Ikebe et al. 
(1999) have revealed that the CEE of Centaurus actually 
gets weaker as the energy increases, but it does not van- 
ish even in the highest ASCA band; the relative promi- 
nence of the central excess is roughly the same between 
the 4.5-6.1 keV and 7.1-10 keV bands (figure 4 of Ikebe 
et al. 1999). Therefore, there must exist central excess 
pressure, which is strong enough to cause the hard-band 
CEE against the displacement by the CCC. The stronger 
CEE in lower energies (< 4 keV) can be attributed to 
the additional contribution from the CCC. Xu (1998) 
has derived a very similar results from 3 A 0335-1-096, 
which also host a strong CCC. 

Thus, from more than half a dozen cD clusters with 
soft X-ray CEE, we have detected the CEE also in hard 
X-rays with energies above 3 keV. Since they are typical 
CEE objects, we can generally ascribe the CEE phe- 
nomenon primarily to a central excess pressure. The 
central ICM temperature decrease (or the CCC) is con- 
cluded to contribute partially to the CEE only in soft 
X-rays, and only in limited objects such as Centaurus 
and 3A 0335-1-096. In terms of the 2T picture developed 
in Section 2, we can rephrase that the hot-phase emis- 
sivity profile of these cD clusters exhibit a central excess 
above the /3-model distribution. These novel findings 



drastically renew our understanding of the CEE phe- 
nomenon. For example, the strong CEE observed from 
the core regions of several distant clusters (e.g. Schindler 
et al. 1997), which is usually attributed to CFs, may in 
reality be the potential shape effect, since these objects 
do not necessarily exhibit strong CCC. 

Figure ^ compares the CCC luminosity of repre- 
sentative objects in the 0.5-3 keV band measured with 
ASCA (taken from table 1), against the CEE lumi- 
nosity Lc calculated in the same band based on the 
ROSAT data analysis by Mohr et al. (1999). Thus, 
Lc considerably exceeds Lc, even though the objects 
plotted here are known to host a rather strong CCC. 
This result can be explained in the following way. Gen- 
erally, we may write the overall luminosity from the 
cluster center region in a given energy band AE as 
l'^\aE) + l'^\aE)+Lc{AE), where l'^\AE) refers 
to the cluster-wide hot emission of which the radial pro- 
file is expressed with a /3 model, while L^^\aE) de- 
scribes the central excess in the hot phase. In this three- 
term expression, the sum of the first two terms gives the 
hot-phase luminosity Li^^AE) from the cluster center re- 
gion, while the sum of the last two terms yields Le{AE). 
Then, the implication of figure |] is that L'i^\AE) of 
these objects considerably exceeds Lc{AE) even in the 
soft X-ray band. This reinforces our previous inference 
made in subsection 2.2, that the cluster core volume is 
dominated by the hot phase. In the hard X-ray band 
(e.g., 3-10 keV), we obviously expect Lc to become neg- 
ligible compared to Lq. 

3.2. Excess central mass as the origin of CEE 

To explore the implication of the central excess pres- 
sure, let us remember that the pressure p{R) of hy- 
drostatic ICM confined by gravity satisfies the relation 
(Sarazin 1988) 

kTR f dlnp 



M{R) 



(10) 



Giinip \dlni?^ 

where M{R) is the total gravitating mass within R, G is 
the gravitational constant, and definitions of k, mp, and 
fi are the same as in equation (^). When the pressure 
profile exhibits a central excess and the ICM is approx- 
imately isothermal, we may write 

p{R) ^ Po{R)aR) , (11) 

where po (R) is the pressure distribution described with 
a P model (or its equivalent, having a flat core), and 
^ (R) > 1 is a non-dimensional factor describing the CEE 
which tends to unity for large values of R. Substituting 
equation ( p] ) into equation (^0|), we obtain 

M{R) = MoiR) + SM{R) . (12) 

Here Mq = —{kTR/Giimp){dlnpo/d\nR) denotes the 
"unperturbed" total gravitating mass distribution ex- 
pressed with a King-type solution having a fiat core, 
while 

kTR /dln^ 



SM{R) = - 



G/XTTip \dlni? 



(13) 
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is an excess mass associated with the central excess pres- 
sure. 

NumericaUy, the ROSAT (Allen, Fabian 1994) and 
ASCA (Ikebe et al. 1999) results both indicate an 
approximate radial dependence of n{R) oc R^^, over 
the R = (10 - 100) h^^ kpc region of the Centau- 
rus cluster. This in turn implies p{R) oc R~^ consid- 
ering the dominance of the hot component. In con- 
trast, the outer-region potential of Centaurus can be 
described by a King-type profile with a core radius of 
i?^ - 7' = 140 h^^ kpc (Matilsky et al. 1985; Ikebe 
et al. 1999), inside which p{R) would be constant 
if there were no CEE. From these facts we may ap- 
proximate as £,{R) oc R^^ in equation (p^), to obtain 
— dln^(i?)/dlni? ^ 1. Taking this as a representative 
case, we can numerically write as 

A more exact treatment considering the cool phase is 
found in Ikebe et al. (1999). 

We can thus infer that a cluster with a hard X-ray 
CEE hosts an excess mass SM of order lO^'^ Mq at 
its center, above the King-type mass distribution. This 
6M produces an additional potential drop at the center, 
which is superposed upon a flat core of the King-type 
potential profile. The potential drop in turn attracts an 
excess amount of hot ICM, thus producing the central 
excess pressure, and hence the hard X-ray CEE (Mak- 
ishima 1995, 1998; Ikebe et al. 1996, 1997b; Xu et al. 
1998). 

The detections of CEE in hard X-rays are a second 
reason why the ICM mass deposition rate was previous 
overestimated (subsection 2.5). That is, the deprojec- 
tion analysis usually proceeds by calculating the ICM 
density profile directly from the observed X-ray surface 
brightness profiles via equation (^, and when the spec- 
tral information is inadequate, the temperature profile 
is derived assuming an ICM pressure equilibrium in a 
plausible gravitating mass distribution, by using equa- 
tion (p^. Therefore, if the central excess mass (or po- 
tential drop) is not properly taken into account or under- 
estimated, the CEE immediately leads to a false central 
temperature decrease, and hence to an overestimate of 
M. 

3. 3. Hierarchical gravitational potential 

The excess mass 5M thus found at the center of a 
cD cluster can most naturally be attributed to the total 
gravitating mass associated with the cD galaxy (Mak- 
ishima 1998), because the CEE is nearly always centered 
on the position of the cD galaxy (Lazzati, Chincarini 
1998) even if it is somewhat offset from the dynamical 
cluster center, and because the excess mass described 
with equation (|lj) is reasonable as the total mass of a 
giant elliptical galaxy. In other words, the gravitational 
potential exhibits a hierarchy^ between a wider compo- 
nent associated with the entire cluster, and a narrower 
one associated with the central galaxy. Such a picture 



has already been considered by various authors (e.g., 
Thomas et al. 1987) in their deprojection analysis, but 
has not been recognized generally as the origin of the 
CEE. 

The cluster/galaxy hierarchy in the total mass distri- 
bution has first been recognized clearly in ASCA data 
of the Fornax cluster and its cD galaxy NGC 1399 by 
Makishima (1995, 1996), Ikebe (1995), and Ikebe et al. 
(1996). These authors successfully reproduced the radial 
X-ray brightness profile observed with ASCA, employing 
a "double-/?" emissivity distribution of the form 

e{R- AE) = eoiR- AE) + ei{R; /:iE) , (15) 

where eg and Ci, both employing three-dimensional (3 
models, describe the cluster component and the cD 
galaxy component, respectively. This decomposition is 
essentially the same as what we have done in subsection 

3.i(L,=Lr-HL«). 

In practice, we fit the observed radial X-ray surface 
brightness profile by a sum of two 2-dimensional (3- 
model components, which are obtained by analytically 
projecting equation (|l5|) onto two dimension. This fixes 
eo(-R) and ei(i?), because a two-dimensional (3 model 
is uniquely related to its three-dimensional counter- 
part. Once the emissivity profile of equation (^5|) is 
thus determined, we can convert it into density profile 
via equation (^, and further into pressure profile as- 
suming approximate isothermality. Obviously, the fac- 
tor ^ in e quation (|ll|) can be given in this case as 
^{R) = + ei{R)/€o{R), and equation (^ shows 
that £1 is responsible for the central excess mass. From 
the derived pressure profile, we can calculate the gravi- 
tating mass profile M{R) via equation (|l0|). 

The mass profile of the Fornax cluster, thus calcu- 
lated by Ikebe et al. (1996), is reproduced in figure ^a. 
There, the 3-dimensional radius is normalized to "inter- 
face radius", i?iF = 72 ± 14 kpc, which is defined as the 
cross-over point of the two emissivity terms in equation 
(|l5|). The mass curve reveals a "shoulder" like struc- 
ture at i? ~ Ryf- Outside Ryp, we observe the cluster- 
scale mass distribution, corresponding to the first term 
of equation (|l^). Inside R ^ Rif, in contrast, there ap- 
pears an additional mass component due to the second 
term of equation (p^, which can be identified with the 
excess mass associated with the cD galaxy. Thus, Rif 
may be regarded as an interface between territories of 
the cD galaxy and the overall cluster (Makishima 1996). 
Substitution of i? = Rif 72 kpc and T = 1.1 keV 
into equation (|l^ yields SM ~ 3 x 10^^ Mq, in rough 
agreement with the data point in figure |]]b. 

The same double-/? modeling has been applied suc- 
cessfully to the ASCA CIS radial profiles of the Hydra-A 
cluster (Ikebe et al. 1997b), Abell 1795 (Xu et al. 1998), 
Abell 2199, AWM7, and 3A 0335-f096 (Xu 1998), aU 
exhibiting the clear hard X-ray CEE (subsection 3.1). 
For several of these objects, we show in figure |5]a the 
rescaled mass curves derived in the same way as for the 
Fornax cluster. The results on the Centaurus cluster 
(Ikebe et al. 1999) properly takes into account the ef- 
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feet of CCC by modifying equation ([T^). Thus, the mass 
curves again bear the characteristic feature at i? ^ ^if- 
Although in some cases (e.g., the Centaurus cluster; see 
subsection 3.5) the apparent slope change in the mass 
curve at i? i?iF might be an artifact caused by the 
particular form of equation (|l5|), the feature is real in 
Abell 1795 (Xu et al. 1998) as weU as in Fornax. There- 
fore, the hierarchical potential shape is suggested to be 
ubiquitous among these objects. 

The potential hierarchy is observed from poorer sys- 
tems as well. The best example is the giant elliptical 
galaxy NGC 4636, around which a large-scale (up to 
~ 300 kpc) ambient X-ray emission was detected with 
ROSAT (Trinchieri et al. 1994). Using ASCA, Mat- 
sushita (1997) and Matsushita et al. (1998) have shown 
that the large-scale emission in fact comes from hot gas 
trapped in the large-scale potential of an "optically dark 
group" , for which NGC 4636 plays the role of "mini-cD" 
galaxy. The X-ray surface brightness was described suc- 
cessfully by equation (jl^), and the calculated mass curve 
in figure ^a again bears the same feature at Rip ~ 23 
kpc, which is not an artifact caused by the particular 
model form of equation ( p^ ) (Matsushita et al. 1998). 
There is evidence (Matsushita 1997; Matsushita et al. 
2000) that X-ray luminous elliptical galaxies can gener- 
ally be regarded as mini-cD galaxies of some larger-scale 
potential structures, whereas X-ray dim ellipticals lack 
such outer potential envelopes. The "X-ray overlumi- 
nous elliptical galaxies" (Vikhlinin et al. 1999) may be 
similar in concept to our "mini-cD" galaxies. 

Figure summarizes the values of Rip and the total 
gravitating mass within it, which are used to normal- 
ize the mass curves in figure [s^a. The results on For- 
nax (Ikebe et al. 1996) and NGC 4636 (Matsushita et 
al. 1998) assume the distances of 20 Mpc and 17 Mpc, 
respectively, which are consistent with the Hubble con- 
stant of 75 km s~^ Mpc~^ {h^o = 1-5). For consistency, 
in figure ^ we therefore employ /150 = 1.5 for the re- 
maining three more distant objects, unlike elsewhere in 
this paper. We see a tight positive correlation between 
i?iF and the enclosed mass. 

The nested two-component emissivity profiles have 
been found even in the ROSAT data. Mulchaey, Zablud- 
off (1998) report that such an X-ray emissivity profile is 
common among galaxy groups. Mohr et al. (1999) suc- 
cessfully employed the double-/? modeling to express ra- 
dial X-ray profiles of a large number of clusters observed 
with ROSAT. Furthermore, the latest Chandra data of 
the Hydra-A cluster clearly reveal a hierarchical surface 
brightness distribution with Rip ^ 100 kpc (David et 
al. 2001). These reports suggest that the hierarchi- 
cal potential shape is ubiquitous among self-gravitating 
systems, although the effects found in the ROSAT data 
may be partially attributable to the ICM temperature 
gradient. 

We can strengthen our discovery of the hierarchical 
potential structure around cD galaxies by referring to 
several optical results suggestive of the same effect. One 
is the large extended stellar envelope of cD galaxies, of 
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which the surface brightness profile deviates significantly 
from the standard de Vaucouleurs law (e.g. Schombert 
1986; Johnstone et al. 1991); the stellar envelope sug- 
gests the presence of a larger-scale potential structure 
surrounding the cD's own potential. The other is the 
outward increase in the stellar velocity dispersion of a cD 
galaxy, which led Dressier (1979) to model the gravita- 
tional potential of a cD cluster as a sum of hierarchically- 
nested three King models. Our hierarchical potential 
structure is essentially identical to the Dressier 's mod- 
eling, although we need only two spatial components 
instead of three. 

As shown so far, the CEE often (if not always) results 
from a nested concentric hierarchy in the gravitational 
potential, formed by the cluster and its cD galaxy. We 
further discuss its implication in subsection 3.5. 

3.4- Central cusp in the gravitational potential 

The CEE is known to be very weak or nearly absent 
in some clusters (Jones, Forman 1984), e.g. the Coma 
cluster, Abell 400, and Abell 1060, which mostly lack 
cD galaxies. Then, what is the gravitating mass dis- 
tribution in these non-cD clusters? Does a King-type 
approximation give an adequate description? 

This issue has been investigated by Tamura et al. 
(1996; 2000) and Tamura (1998), through detailed X- 
ray studies of Abell 1060. Being a prototypical weak- 
CCC object with a nearly isothermal ICM as shown in 
subsection 2.3, this cluster is also known to have lit- 
tle CEE (Jones, Forman 1984). However, the ROSAT 
PSPC radial brightness profile exhibits a weak yet sig- 
nificant excess above a /3-model (Tamura 1998; Tamura 
et al. 2000), and hence cannot be explained by emis- 
sion from an isothermal ICM confined in a King-type 
potential. 

In order to reconcile the good isothermality indicated 
by the ASCA data and the weak CEE found in the 
ROSAT PSPC data, Tamura (1998) and Tamura et al. 
(2000) have resorted to employ, instead of a King-type 
model, a total mass density distribution of the form 

where pg, Rs, and ^ are positive parameters. This 
distribution, found by A'^-body simulations (Navarro et 
al. 1996; Fukushige, Makino 1997; Moore et al. 1998; 
Ghigna et al. 2000), has a singular cusp at the cen- 
ter, but the implied gravitational potential remains finite 
there as long as C < 2. For C = 1) this formula reduces 
to the "universal halo" profile of Navarro et al. (1996). 
Then, the overall ASCA and ROSAT data of AbeU 1060 
have been reproduced successfully in terms of emission 
from a nearly isothermal ICM confined in the potential 
corresponding to equation ( |l6| ) with ( — 1.5± 0.1, while 
the King-model profile has been ruled out (for detail, 
see Tamura et al. 2000). Thus, the X-ray data of this 
prototypical non-cD cluster suggest the gravitational po- 
tential with a central cusp, and the derived value of ( 
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is in a very good agreement with that predicted by the 
iV-body simulation (Ghigna et al. 2000). 

Similarly, Markevitch et al. (1999) reported that the 
mass profiles of Abell 2199 and Abell 496 can be ex- 
pressed well by equation (^6|). In addition, by analyzing 
a large number of objects observed with ASCA, Sato et 
al. (2000) concluded that equation (|l6|) with C = 1 can 
successfully reproduce their dark halo shapes: however, 
details are not given in their paper. 

3. 5. The two types of the gravitational potential shape 
How does the cuspy potential profile, indicated by A^- 
body simulations and observed from Abell 1060 (sub- 
section 3.4), relate to the hierarchically nested halo-in- 
halo type potential profile (subsection 3.3) found among 
cD clusters? Observationally, the two phenomena are 
qualitatively similar, because they both imply an ex- 
cess gravitating mass at the cluster center as noted in 
subsection 3.2. Actually, the hard X-ray GEE of the 
Centaurus cluster can be explained by either of the two 
potential models (Ikebe et al. 1999), because the strong 
CCC makes the potential shape ambiguous. The case of 
Abell 2199 may be similar, as its potential shape can be 
expressed by either the double-/? model (Xu 1998) or the 
cuspy model (Markevitch et al. 1999). Conversely, the 
weak CEE of Abell 1060 may be described alternatively 
by a double-/? model. 

From physics viewpoint, comparison of the central ex- 
cess mass found in subsections 3.2 and 3.3 against the 
cD's stellar mass yields a baryon faction of 0.2 — 0.3, at 
least in Fornax (Ikebe et al. 1996), Centaurus (Ikebe 
et al. 1999), and NGC 4636 (Matsushita et al. 1998). 
Therefore, the central excess mass is still dark, imply- 
ing that the dark matter is clustered on the two distinct 
spatial scales. This concept is essentially the same as 
the physics behind the formation of the cuspy potential, 
that smaller-scale dark matter halos formed in earlier 
epochs survive the subsequent hierarchical merging, and 
remain as subhalos near the cluster center (Navarro et 
al. 1996; Dubinsky 1998; Ghigna et al. 2000). 

In contrast to the above arguments for the close re- 
semblance between the cuspy and hierarchical poten- 
tial models, some observational facts argue against their 
identification. First of all, the ASCA data of the For- 
nax cluster reveals the two distinct spatial scales with- 
out ambiguity. Similarly, the hierarchical structure seen 
in NGC 4636 can be derived from the surface bright- 
ness profile in an non-parametric manner (Matsushita 
et al. 1998), using Monte-Carlo simulations instead of 
the double-/? modeling. Furthermore, the mass curve of 
Abell 1795 is better interpreted as exhibiting a hierar- 
chical structure than a cuspy potential. According to 
Williams et al. (1999), the cluster core mass estimated 
through gravitational lensing tends to exceed those pre- 
dicted by the A^-body simulations, presumably due to 
the excess mass associated with the cD galaxy. We there- 
fore regard the two potential structures as distinct. In 
fact, as illustrated in figure the hierarchical potential 
can produce a much stronger CEE than the cuspy po- 



tential can (Makino et al. 1998; Tamura et al. 2000). 
These implications are independent of our 2T formalism. 
The difference between the two potential shapes may be 
consistent with the statistical result by Fujita, Takahara 
(1999), that clusters form a two-parameter family, with 
the two parameters being the system mass and the de- 
gree of central mass concentration. 

One likely scenario is that the central cusp is formed 
essentially in all clusters by the dark matter distribu- 
tion, whereas the cusp develops into a much deeper cen- 
tral dimple as the cD galaxy builds up through, presum- 
ably, under a considerable baryonic energy dissipation, 
as evidenced by the factor 2-3 higher baryon fraction 
in the central dimple region than in the entire cluster. 
To assess such a possibility, in figure [s^ b we draw a 
straight line which indicates the scaling relation found 
by Navarro et al. (1996) between i^max and the total 
integrated mass contained within it, where i?max is the 
radius at which the circular velocity of the halo becomes 
maximum. Although the observed data points lie by a 
factor 2-3 above the prediction, they will lie closer to 
the line, if we subtract the mass contribution from the 
underlying cluster component. Therefore, the central 
excess mass itself may be regarded as a self-gravitating 
dark halo, and a cD cluster can be regarded as a halo- 
in-halo system. We expect that future numerical simu- 
lations, fully taking into gas dynamics, will be able to 
reproduce the hierarchical potential profile found in cD 
clusters. 

In any case, our measurements of the cluster poten- 
tial profile imply the presence of a particular struc- 
ture in the dark matter distribution underlying each cD 
galaxy. This disagrees with the scenario (e.g., Fabian 
1994) which describes cD galaxies as a result of gas con- 
densation in CFs. 

4. ICM METALLICITY IN THE CLUSTER 
CENTER REGIONS 

Production, confinement, and transport of heavy ele- 
ments are a third important aspect of the X-ray study 
of clusters of galaxies. ASCA observations have for the 
first time enabled systematic studies of spatial distribu- 
tions of heavy elements in the ICM, particularly iron and 
silicon. We again find characteristic phenomena around 
cD galaxies. 

4.I. Cluster-wide properties 

Excluding the CCC regions, the average iron abun- 
dance of ICM is 0.2-0.3 solar, without significantly de- 
pending on the cluster richness (Fukazawa et al. 1998). 
As a result, the total iron mass in the ICM of a cluster 
amounts to (1 — 3) x 10"'^ of the total stellar mass therein 
(Tsuru 1992; Arnaud et al. 1992). Equivalently, the iron 
mass to light ratio (IMLR; Ciotti et al. 1991; Renzini et 
al. 1993; Renzini 1997), i.e. the iron mass in the ICM 
normalized to the stellar light, becomes (0.5 — 3) x 10~^ 
in the solar unit. This amount of iron is indeed com- 
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parable to the total iron locked in the stellar interior 
(Renzini et al. 1993). 

Such a large amount of iron (and other heavy ele- 
ments) in the ICM must have been produced in the 
stellar interior of member galaxies, particularly ellipti- 
cals (Arimoto, Yoshii 1987), and subsequently ejected 
into the intra-cluster space. However, observations with 
Ginga (Awaki et al. 1991), ROSAT (Forman et al. 
1993), BBXRT (Serlemitsos et al. 1993), and ASCA 
(Awaki et al. 1994; Loewenstein et al. 1994; Matsumoto 
et al. 1997) have failed to detect the expected trace of 
such metal ejection processes around individual galax- 
ies: the ISM metallicity of X-ray luminous ellipticals is 
at most 1 solar, and those of X-ray fainter ones can be 
even lower (Matsushita et al. 1997, 2000). While these 
values are consistent with those expected from stellar 
mass-loss, there is not much room left for the metal en- 
richment of the ICM by supernova (SN) products. This 
clearly indicates that the SN products have been re- 
moved quickly from individual galaxies, and transported 
into the intra-cluster space, via, e.g., ram-pressure strip- 
ping and energetic outflow. Indeed, the iron mass now 
contained in the ISM of all member galaxies of a clus- 
ter would sum up to make only a few percent of that 
contained in the ICM. 

In order to reinforce this view, in figure |7] we summa- 
rize the IMLR of various objects, as a function of their 
plasma temperature serving as a measure of the sys- 
tem richness. The figure, originally devised by Fukazawa 
et al. (1996) and Ishimaru (1996, 1998), includes clus- 
ters with different richness, galaxy groups, and elliptical 
galaxies. Thus, the IMLR clearly decreases as the sys- 
tem gets poorer. Evidently, poorer systems have lost 
most of the heavy elements produced in them, presum- 
ably because they have too low an efficiency of grav- 
itational confinement of the metal-enriched SN prod- 
ucts (Fukazawa 1997; Matsushita 1997; Fukazawa et al. 
2000). The heavy elements are thus inferred to be escap- 
ing extensively from objects of lower hierarchy, and the 
escaped materials enrich the systems of higher hierarchy. 

In figure |7], elliptical galaxies form two distinct sub- 
groups, as discovered by Matsushita (1997) and Mat- 
sushita et al. (2000). One class ("X-ray extended" ones) 
comprises X-ray luminous objects such as NGC 4636, 
which also have higher ISM abundances. They exhibit 
extended X-ray morphology indicative of larger-scale 
outer potential envelopes, for which they act as mini- 
cD galaxies (subsection 3.3). The other class ("X-ray 
compact" ones) consists of ellipticals with low X-ray lu- 
minosities (< 10^° ergs s~^), which lack outer poten- 
tial envelopes. The difference in IMLR between the two 
types can also be understood in the context of metal 
confinement and escape, because the outer potential en- 
velopes associated with objects of the former class are 
though to improve the confinement efficiency of the en- 
ergetic SN materials (Matsushita 1997). 

The ASCA observations have for the first time yielded 
systematic measurements of the silicon-to-iron abun- 
dance ratio in the ICM, which is a key to determining 
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the relative importance of type la and type II super- 
novae. In the outer regions of medium-richness clus- 
ters, silicon has been found to be systematically over- 
abundant relative to iron by a factor of 1.5-2 in solar 
units (e.g., Mushotzky et al. 1996; Tamura et al. 1996; 
Fukazawa et al. 1998). This reveals the dominant role of 
type II SNe in the process of metal enrichment of ICM. 
However, Fukazawa (1997) and Fukazawa et al. (1998) 
have discovered that the silicon over-abundance gradu- 
ally disappears towards poorer clusters, and the poorest 
ones exhibit solar- like Si/Fe ratios implying a significant 
contribution from type la SNe as well. In other words, 
"silicon mass to light ratio" of the ICM decreases con- 
siderably, e.g. by a factor of 2, from rich clusters to 
the poorest ones, whereas the IMLR is approximately 
constant among clusters (figure As pointed out by 
Fukazawa et al. (1998), one possible explanation of this 
effect is a selective escape of silicon; presumably, silicon 
has been supplied early in the cluster evolution mainly 
in the form of galactic winds created by type II SNe, 
which were energetic enough to escape from poorer sys- 
tems. In contrast, iron is likely to have been supplied by 
more prolonged activity of type la SNe without forming 
such energetic outflows. 

4-2. Large-scale radial behavior of IMLR 

Radial changes in the IMLR are expected to provide 
further clues as to the metal production, conflnement, 
and escape. For this purpose, we extend the concept of 
IMLR into "radial IMLR profile" , defined as 

$(i?) = A/Fc(i?)/i*(i?) , (17) 

where M-pc{R) and L^{R) denote the radial profile of 
the iron mass in the ICM and that of the stellar light, 
respectively, both integrated within R and expressed in 
solar units. In outer regions of the cluster, L^{R) gen- 
erally increases more slowly with R than the integrated 
ICM mass profile, Micm(^)- Taking the AWM7 clus- 
ter for example, the differential forms of its Micm(^) 
and L^{R) can be approximated by /? models of /3 = 
0.58t^;^^ and /? - 0.8, respectively (Ezawa et al. 1997). 
Therefore, at large radii, the ratio Micm{R) / L*{R) in- 
creases as oc where 0.66 = 3 x (0.8 - 0.58) 
is the logarithmic slope difference between Micm(-R) 
and L*(i?) in outer regions. We then expect <&(i?) oc 
i?0-66z(i?), where Z{R) is the radial profile of the ICM 
iron abundance. 

In many clusters, Z(R) has so far been considered 
approximately constant as a function of R, except in 
the central regions. However, Ezawa et al. (1997) have 
discovered that Z{R) of AWM7 in fact decreases over 
a cluster-wide spatial scale as Z{R) oc out- 
side i? 4' = 120 h'^Q kpc; Mpe is less extended than 
MicM- This effect just cancels out the outward increase 
in Micm{R)/L^{R), and makes $(-R) almost radially 
constant at 2.8 x IQ-^. Namely, Mp^iR) of AWM7 be- 
haves approximately proportional to (R) over spatial 
scales from ~ 100 kpc to ~ 1 Mpc, implying that the 
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iron density in the ICM faithfully traces the stellar light 
distribution. 

Similar results have been obtained with ASCA from 
the Perseus cluster (Ezawa 1998), Abell 4059 (Kikuchi e 
al. 1999), and Abell 2029 (Molendi, De Grandi 1999) as 
well. We presume that these effects, i.e., the large-scale 
outward decrease in Z{R) and the associated constancy 
of $(i?) outside R ^ 100 kpc, are generally present in 
many clusters. The reason why the effects have so far 
been observed from the four particular clusters may be 
ascribed to their high X-ray surface brightness and large 
angular extent. With the future instrumentation, we 
hence expect this effect to be detected ubiquitously from 
many clusters. Actually, the XMM-Newton observation 
of the southern cluster Sersic 159-03 clearly reveals a 
large-scale ICM metallicity gradient over a radial scale 
of - 170 kpc (Kaastra et al. 2001). 

These results provide one of the first direct confirma- 
tions of the general consensus that the heavy elements 
in the ICM were indeed ejected by the cluster mem- 
ber galaxies into the intra-cluster space. Furthermore, 
the constant IMLR profiles suggest that the metals do 
not travel over a large distance (Ezawa et al. 1997), 
even though they must be removed quickly from their 
source galaxies (subsection 4.1), via, e.g., ram pressure 
stripping and energetic outfiow. Presumably, the galax- 
ies have been swimming in the ICM while continuously 
depositing metal-enriched gas (e.g., Charlton, Salpeter 
1989). The metals then "mixed" into the ICM and soon 
became hydrostatic, resulting in a radial distribution 
which is similar to the galaxy distribution in a statis- 
tical sense. 

We must remark here that the observations do not 
necessarily exclude homologous expansion or contraction 
in Micm(^) or i*(i?), in such a way that their core radii 
change but their /3's remain unchanged. In these cases, 
^(i?) would vary over the core region, but would still 
remain constant outside it. This urges us to study the 
metallicity behavior in the cluster center (< 100 kpc), 
in the following two subsections. 

4-3. ICM metallicity near the cluster center 

The ICM metallicity has been observed to increase sig- 
nificantly toward the center of many cD clusters, includ- 
ing Virgo (Koyama et al. 1991; Matsumoto et al. 1996), 
Centaurus (Fukazawa et al. 1994; Ikebe 1995; Ikebe et 
al. 1999), AWM7(Xuetal. 1997), Abell 496 (Hatsukade 
et al. 1998), Abell 262 (David et al. 1996), AbeU 4059 
(Kikuchi et al. 1999), and possibly several more objects 
(Ohashietal. 1994). We can also add recent results from 
Chandra on the Hydra- A cluster (David et al. 2001) and 
XMM-Newton measurements of Abell 1795 (Tamura et 
al. 2001). This effect is limited to the central 100 
kpc scale where the CCC and CEE phenomena are also 
observed, and much more localized to the center than 
the outward decrease in Z{R) described in subsection 
4.2. Although these results are derived assuming the two 
ICM phases to have the same abundances, the metallic- 
ity increase in the hot phase is a robust result, as shown 



by the enhanced Fe-K lines in the central regions of most 
of these objects. Furthermore, the CCC of the Centau- 
rus cluster is at least as metal enriched as the central 
hot phase (Ikebe et al. 1999). Therefore, we presume 
that the central increase in Z{R) in these objects occurs 
in both hot and cool phases. 

In contrast to these cD clusters, Abell 1060, the typi- 
cal non-cD cluster, exhibits an essentially constant Z{R) 
from the center up to ~ 250 h'^Q kpc (Tamura et al. 
1996). The statistical study of 40 clusters by Fukazawa 
(1997), described in subsections 2.2 and 2.3, confirm that 
this can be regarded as a systematic difference between 
cD and non-cD clusters: the central ~ 100 kpc region 
of cD clusters exhibit systematically higher abundances 
than their outer regions, while non-cD clusters exhibit 
spatially uniform ICM abundances over the core region. 
(In the peripheral regions, both types of clusters are ex- 
pected to exhibit the large-scale abundance decrease de- 
scribed in subsection 4.2.) We therefore conclude that 
the ICM abundance generally increases within ^ 100 
kpc of cD galaxies. 

The above conclusion is reinforced by figure where 
we plot the central iron abundance against the Qc/Qh 
ratio after Tamura et al. (1997). This results is based 
on a subsample of table 1 with redshifts < 0.04, because 
the central metal-enriched regions of more distant ob- 
jects are difficult to resolve with ASCA. Thus, the two 
quantities exhibit a tight positive correlation, with Cen- 
taurus and Virgo at one end of the distribution while the 
non-cD clusters (AbeU 1060, AbeU 400, and AbeU 539) 
at the other end. Since the CCC is a clear signature 
of cD galaxies (subsections 2.4 and 2.6), the correlation 
indicates a close relation between the presence of a cD 
galaxy and the central metallicity enhancement. A par- 
ticular advantage of figure |]] is that its implication is 
independent of the definition of cD and non-cD clusters. 

In a typical cD cluster, the excess iron mass contained 
in the central region (hot plus cool phases) amounts to 
- 7x 10^ Mq or less (Fukazawa et al. 2000). This can be 
supplied over the whole lifetime of a single giant ellipti- 
cal galaxy via type la SNe (Renzini et al. 1993). There- 
fore, the excess metals at the cluster center can most 
naturally be regarded as a product of the cD galaxy, as 
argued by several authors (e.g., Fukazawa et al. 1994, 
2000; Ikebe et al. 1999; Kikuchi et al. 1999; David et al. 
2001). The ejected metals are thought to have mostly re- 
mained in the cluster core region, because the cD galaxy 
is free from ram-pressure stripping effects unlike the 
other galaxies that are moving through the ICM. This 
viewpoint is reinforced by the systematic difference seen 
in the ICM chemical composition between central and 
outer regions of cD clusters (Fukazawa 1997; Fukazawa 
et al. 2000). Namely, in the metal-enriched cluster cen- 
ter regions (again assuming the hot and cool phases to 
have the same abundances), the Si/Fe ratio is approxi- 
mately unity in solar units on average, even though the 
silicon over-abundance is observed from outer regions 
particularly when the cluster is relatively rich (subsec- 
tion 4.1). The same conclusion has been derived by 
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David et al. (2001) on the Hydra-A cluster using Chan- 
dra. This solar-like Si/Fe ratio resembles those mea- 
sured in the ISM of X-ray luminous elliptical galaxies 
(Matsushita 1997; Matsushita et al. 1998, 2000), and 
suggests a significant contribution from type la SNe in 
the cD galaxy. This particular result, considered to hold 
both for the hot and cool phases, also gives a convincing 
support to our view presented in subsection 2.6 that the 
CCC is associated with the cD galaxy. 

In contrast to cD clusters, the central regions of non- 
cD clusters are usually populated with several galaxies 
of comparable luminosities. Since none of such galax- 
ies is at rest in the gravitational potential, they are all 
subject to the ram-pressure stripping, and metals they 
produced must have been transported efficiently to the 
intra-cluster space just like those from the other mem- 
ber galaxies. Furthermore, the relative motion of these 
galaxies will efficiently mix the metal-rich ejecta into 
the ICM. These ideas, already pointed out by Tamura 
et al. (1996), give a natural account of the spatially 
uniform abundances of non-cD clusters. This mecha- 
nism may apply even to some cD clusters. For example, 
Tamura et al. (1996) argue that Z{R) of the Fornax 
cluster stays rather constant at the center, because its 
cD galaxy, NGC 1399, has a close companion NGC 1404 
which must be in a strong interaction with NGC 1399. 

Using ROSAT images and spatially averaged ASCA 
spectra, Allen, Fabian (1998) showed that the emission- 
weighted ICM abundances of "cooling-flow (CF)" clus- 
ters are systematically higher than those of "non cooling- 
flow (NCF)" objects. The sample they used (rather 
distant objects) is mostly disjoint from ours, and their 
subsample classification is different from ours (cD vs. 
non-cD). Nevertheless, among their 21 "CF" clusters, 
five have Bautz-Morgan (B-M) type I, one has type I- 
II, and four have type II; none is classified as B-M type 
II-III or III, while the B-M type is unavailable for 10 
objects. On the contrary, of their nine "NCF" clusters, 
one has B-M type II, three have type II-III, four have 
type III, and one unknown, but none is type I or I-II. 
Therefore, the Allen & Fabian's classification is nearly 
identical to ours. Their "CF" clusters in fact contain 
cD galaxies, and the associated central metallicity en- 
hancements, augmented by the GEE, presumably made 
the emission-weighted metallicity higher than those of 
"NCF" ones. Thus, their results are consistent with 
ours, but we consider the ICM metallicity to be more 
closely related to the presence/absence of a single domi- 
nant galaxy at the cluster center, rather than to the CF 
strength. 

In an attempt to reconcile the abundance increase of 
the Centaurus cluster with the CF hypothesis, Reiseneg- 
ger et al. (1996) argued that the metals ejected from 
the cD galaxy are swept back by CF to become com- 
pressed at the center. They hence predict an anti- 
correlation between the central increment in Z{R) and 
the CF strength, because the metals would drop out of 
the ICM as the cooling proceeds. However, what we 
observe in figure ^ is quite opposite to the prediction. 



and rules out their interpretation. Therefore, radiative 
metal dropouts are unlikely to be playing a major role. 

4- 4- IMLR profiles at cluster centers 

How does ^(i?) behave near the center of clusters? For 
this purpose, in figure |] we show $(i?) of Abell 1060, 
calculated from Tamura et al. (2000). Thus, $(i?) 
clearly decreases to the center. This is a trivial con- 
sequence of the two well established facts, that the stel- 
lar light profile is more centrally peaked than the ICM 
mass profile, and that the ICM metallicity of Abell 1060 
is spatially uniform (subsection 4.3). 

The central decrease in $(i?) might be a result of cool- 
ing dropouts of the metals. To examine this possibility, 
in figure ^ we have plotted <&(i?) of the Centaurus clus- 
ter as well, taken from Ikebe et al. (1999). Although 
we again observe a central drop in the IMLR profile, it 
is milder than that of Abell 1060, because of the strong 
metallicity increase at the center. Thus, the Centaurus 
cluster that has a much higher CF rate exhibits a less 
marked decrease in $(i?) than Abell 1060 that has little 
CF. Therefore, the radiative cooling effects cannot be 
the main cause of the decrease in 3>(i?), in agreement 
with our inference made in subsection 4.3. Incidentally, 
the curves of ^(i?) in figure ^ are not affected signif- 
icantly by the insufficient angular resolution, because 
they have both been calculated by fully taking into ac- 
count the instrumental point spread function. 

The central decrease in IMLR is observed from other 
clusters as well. For example, at a representative radius 
of i? - 100 kpc, AWM7 and Abell 4059 exhibit 
IMLR of 3.2 X 10-3 (converted from Xu et al. 1997) 
and 4 x 10"^ (Kikuchi et al. 1999), respectively. Both 
these values are very close to those of Abell 1060 and the 
Centaurus cluster. From this result and those obtained 
in subsection 4.2, we can generally characterize $(i?) by 
two features, namely its constancy outside 100 200 
kpc, and its decrease toward the center. We may restate 
the idea that L^{R) and M-pe{R) share nearly the same 
value of (3 which in turn is larger (steeper) than that for 
-Wicm(-R), while L^{R) has a smaller core radius than 

AfFc(i?). 

Because the ejected metals are likely to be conserved 
in the ICM, the suggested difference in the core radius, 
between ^^(i?) and MFe(-R), is thought to reflect time- 
dependent changes in the spatial distributions of the 
stars and/or metals. More specifically, the phenomenon 
can be explained as a result of either a radial contrac- 
tion in the stellar distribution, or a radial expansion in 
the metal distribution. The latter mechanism must be 
operating to a certain extent, because we have invoked 
repeatedly the metal escape effects. However, the uni- 
form ICM abundances at the center of non-cD clusters 
would require that the suggested metal outflow is bal- 
anced in detail by a radial expansion of the primordial 
ICM. Since such a fine tuning is generally unlikely, we 
consider the former mechanism, i.e. the galaxy infall, 
to be more dominant. The galaxies, while continuously 
ejecting metals, are therefore suggested to have gradu- 
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ally fallen to the cluster center, or merged into bigger 
ones residing in central regions, to achieve the centrally- 
peaked distribution of the stellar component. 

The values of IMLR so far attributed to the cluster 
core regions, several times 10~^, all refer to the iron 
contained in both hot and cool phases. In contrast, if 
we consider the CCC only, the IMLR further reduces to 
several times 10"'*. This is comparable to those of the 
ISM of X-ray luminous ellipticals. This gives further 
support to our conclusion made in subsection 2.6 that 
the CCC is the cD's ISM. Presumably, a major fraction 
of the iron produced in the cD galaxy is contained in the 
hot phase in the form of a central metal excess, which 
represents what would have normally been removed if 
the galaxy were moving through the ICM. 

5. DISCUSSION 

5.1. Summary of new results 

So far, many cluster X-ray investigators purposely 

avoided analyzing the data from cluster central regions, 
except when they try to develop the CF scenario. We 
have challenged this issue through ASCA observations 
of nearby clusters, and arrived at a new picture of the 
gas and mass structure of the central regions of clus- 
ters. Our results are summarized in the following three 
aspects. 

In section 3, we have found that the gravitational po- 
tential exhibits a marked deepening at the cluster cen- 
ter, as evidenced by the presence of CEE in hard X-rays. 
More specifically, the gravitating mass distribution can 
be described by either a hierarchical model involving two 
characteristic scale lengths (corresponding to the cluster 
and the cD galaxy), or a scale- free distribution with a 
central cusp as predicted by numerical simulations. We 
tentatively assign the former type of distribution to cD 
clusters, and the latter to non-cD ones, although the 
sample is still small. In particular, we emphasize the 
halo-in-halo structure formed around each cD galaxy. 

The study in section 2 employing the 2T formalism 
has lead us to propose that the CCC, which is specific to 
cD clusters, has characteristics of the ISM (inter-stcllar 
medium) of the cD galaxy filling the cD's own potential, 
rather than of a cooling portion of the ICM. To the five 
supporting facts listed in subsection 2.6, we can add the 
two metallicity arguments presented in subsection 4.3, 
that the regions around cD galaxies exhibit metallicity 
enhancements, and that the chemical compositions there 
somewhat differ from those of the bulk ICM. The effect 
of ICM cooling, if any, turns out to be much smaller 
than was thought previously, suggesting that significant 
ICM heating mechanisms are operating. 

Our study of the ICM metallicity in section 4 incor- 
porating the concept of the IMLR profile have revealed 
two interesting inferences. One is that the metals have 
been escaping extensively from a poor system, a part of 
which can be trapped by the potential of a surrounding 
system of a higher hierarchy. The other, though more 
speculative, is that the radial galaxy distribution has 



been gradually shrinking relative to the ICM distribu- 
tion, and possibly relative to the dark-matter distribu- 
tion as well. The observed decrease in the IMLR profile 
toward the center may be explained by a combination of 
these two efi^ects. 

As a straightforward application of our results, we can 
for the first time solve the long-lasting confusion as to 
the nature of the bright X-ray emission associated with 
M87, the cD galaxy of the Virgo cluster. So far, some 
authors (e.g., Takano et al. 1989; Bohringer et al. 1994) 
interpreted it as a part of the Virgo ICM emission, while 
others (e.g., Fabricant, Gorenstein 1983; Beuing et al. 
1999) discussed it in terms of the ISM emission from 
M87 as an elliptical galaxy. According to our picture, 
the central cool phase corresponds to the ISM associated 
with M87, whereas the hot phase around M87 represents 
the Virgo ICM of which the density is enhanced by the 
potential drop associated with M87. This statement, 
first suggested by Matsumoto et al. (1996), concisely 
summarizes the present work. 

A large amount of new data are accumulating from 
Chandra and XMM-Newton observations. Although de- 
tailed comparison of these results with ours is yet to be 
carried out, the main points we have emphasized in this 
paper are mostly being reconfirmed and reinforced by 
these two powerful missions. 

5.2. Heating and thermal stability of the cool compo- 
nent 

As mentioned in subsection 2.6, our scenario is subject 
to several theoretical issues yet to be solved; (1) how 
the cool ISM phase is separated and thermally insulated 
from the hot ICM phase, against rapid heat conduction; 
(2) what supplies the cool phase with the large amount 
of energy (up to ~ 10^^ ergs s~^) necessary to sustain the 
X-ray radiation and prevent it from the thermal collapse; 
and (3) how the heating balances the cooling in a stable 
manner. Referring to Makishima (1997ab, 1999ab) and 
Ikcbe ct al. (1999), we below present some speculative 
ideas that might answer these questions. 

The ICM is generally magnetized up to a few micro 
Gauss level (e.g., Kronberg 1994; Eilek 1999). The mag- 
netic field is known to be particularly strong near cD 
galaxies, reaching 10 — 100 /zG (Taylor et al. 1999). As 
a natural consequence, we may invoke magnetic fields as 
the required thermal insulator between the two phases 
(Makishima 1994b, 1997b, 1999a). Even though the 
actual field strengths may fall below the equipartition 
value (~ 30 /K G in typical cases), the heat conduction 
can be suppressed by many orders of magnitude. This 
could solve the first question. 

As to the second question, various heating source 
candidates have long been considered, including super- 
novae, active galactic nuclei, drag due to the galaxy mo- 
tion, and so on. A general consensus is that none of 
these candidates provide heating luminosity high enough 
to sustain the radiative energy loss of the cool phase 
(e.g., Bregman, David 1989; Fabian 1994). However, 
one important possibility has not been fully consid- 
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ered; i.e., magnetohydrodynamic (MHD) effects as the 
member galaxies move through the ICM (Makishima 
1997ab, 1999b). Because the ICM is such an ideal classi- 
cal plasma with extremely high magnetic Reynolds num- 
bers, we expect the galaxy motion to cause significant 
MHD turbulence and frequent magnetic reconnection 
(Norman, Meiksin 1996). Then, the kinetic energy of 
galaxies may be dissipated with a much higher efficiency 
than would be expected when the gas is neutral. The 
recently reported X-ray "wakes" (Drake et al. 2000), 
based on the ROSAT HRI image of Abell 160, may pro- 
vide evidence for such effects. Although calculation of 
the heating luminosity of the proposed mechanism is 
subject to large uncertainties, a crude estimate suggests 
that it might work (Makishima 1999b). 

If the proposed mechanism invoking the MHD effects 
is actually operating, we should expect reactions from 
the ICM to the member galaxies. There are in fact sev- 
eral hints of such effects in the existing data. One is the 
suggestion made in subsection 4.4, that the radial galaxy 
distribution may have been shrinking relative to that of 
the ICM. It might be that the galaxies have lost some 
portion of their kinetic energies through interaction with 
the ICM, and have gradually fallen inwards, while the 
ICM expanded by receiving the released energy. The 
observed increase in the fraction of elliptical galaxies to- 
ward the cluster center (e.g., Whitmore, Gilmore 1993) 
can be another piece of evidence, because the drag force 
will accelerate mergers of smaller spirals into larger ellip- 
ticals. Similarly, the observed galaxy "metamorphosis" 
from distant to nearby clusters (Dressier et al. 1994) 
could be a result of enhanced galaxy interactions and 
mergers through the ICM drag (Makishima 1999a). 

In addition to the proposed galaxy-ICM interaction, 
we can speculate on another ICM heating mechanism lo- 
calized around each cD galaxy, where the heating energy 
is most vitally needed. As described in section 3 and fig- 
ure |]]b, a cD galaxy is likely to form a self-gravitating 
core, immersed in the larger cluster halo. The self- 
gravitating energy reaches ^ G{SM)^ /Rip 1 x 10^^ 
ergs, where SM refers to equation (|lj) and Rif ~ 100 
kpc is the typical interface radius (subsection 3.3). As 
the self- gravitating core gradually shrinks, possibly un- 
der a significant baryonic influence, half the released 
gravitational energy would be radiated, while the re- 
maining half would be spent in the ICM heating. If, 
e.g., ^ 10% of the self-gravitating energy is thus released 
over the Hubble time, the available heating luminosity, 

1 X 10'*'^ ergs s'-^, may be sufficient to prevent the 
ICM from radiative collapse (Makishima 1999b). 

Even though the heating energy may be available, the 
volume cooling rate of the ICM due to radiation is pro- 
portional to as given by equation (^), while the vol- 
ume heating rate is normally proportional to n. There- 
fore, it is usually difficult for the heating mechanisms to 
stably balance the cooling. However, the solar corona, 
confined within many magnetic loops, is thermally stabi- 
lized by the self-regulating Rosner- Tucker- Vaiana (1978) 
mechanism (Kano, Tsuneta 1995). As an analogy, we 



speculate that the magnetic field in the cluster core re- 
gion takes a form of numerous loops anchored to the cD 
galaxy; the loops arc surrounded by the hot phase, and 
their interior is occupied by the cool phase. Makishima 
(1997b) termed this concept "cD corona": such struc- 
tured fields are suggested in the literature, (Ge, Owen 
1994; Taylor et al. 1994; Owen et al. 1999; Bohringer 
1999). If so, the heating and cooling could be stably bal- 
anced by the Rosner- Tucker- Vaiana mechanism. This 
idea might give an answer to the third question. 

6. Conclusion 

By combining several observational works on central 
regions of galaxy clusters achieved with ASCA, we have 
arrived at a novel view therein. This view describes 
the region around a cD galaxy as a site of significant 
and active evolution, where a large amount of heavy ele- 
ments are produced, a self-gravitating core develops, and 
presumably certain ICM heating processes are operat- 
ing. The scenario makes a contrast to the previous view 
which emphasized the role of radiative plasma cooling. 
The recent Chandra and XMM-Newton results clearly 
favor our scenario. 
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Table 1.. Two-temperature modeling of a sample of nearby clusters, quoted from Pukazawa (1997) and Pukaaawa et al. (2000). 
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± 


0, 


.15 


2, 


.04 


± 


0, 


,63 


0, 


,22 


±0.07 


6.3 


2.6 ± 1.4 


A119 


0, 


.044 


II-III 


5, 


.52 


± 


0, 


.27 














< 


0.16 


< 1.2 




A3558 


0, 


.048 


I 


5, 


,17 


± 


0, 


.25 














< 


0.05 


< 1.5 




A2147 


0, 


.036 


III 


4. 


.92 


± 


0, 


.22 














< 


0.20 


< 0.6 




A496 


0, 


.032 


I 


4. 


.12 


± 


0, 


.12 


1, 


.80 


± 


0, 


,14 


0, 


,62 


±0.08 


3.2 


10.5 ± 1.8 


A2199 


0, 


.030 


I 


4, 


,10 


± 


0, 


,07 


1, 


.78 


± 


0, 


,32 


0, 


,17 


±0.05 


1.2 


3.1 ± 1.2 


A4059 


0, 


.048 


I 


4, 


,03 


± 


0, 


,11 


1, 


,38 


± 


0, 


,25 


0, 


,08 


±0.06 


0.5 


0.9 ±0.7 


AWM7 


0, 


.018 


I 


3, 


,76 


± 


0, 


,08 


1, 


,83 


± 


0, 


,46 


0, 


,27 


±0.13 


0.3 


6.0 ±3.5 


Centaurus 


0. 


Oil 


I-II 


3. 


,75 


± 


0, 


.07 


1, 


.37 


± 


0, 


,02 


1, 


,02 


±0.09 


0.8 


12.4 ± 1.2 


MKWSs 


0. 


.043 


II-III 


3. 


,69 


± 


0, 


.20 


1, 


.45 


± 


0, 


,61 


0, 


,09 


±0.07 


0.6 


1.4 ± 1.4 


A2063 


0, 


.034 


II 


3, 


,68 


± 


0, 


,11 


1, 


,64 


± 


0, 


,64 


0, 


,19 


±0.12 


0.4 


3.6 ± 3.0 


A2634 


0, 


.031 


II 


3, 


,64 


± 


0, 


,26 


1, 


.45 


± 


0, 


,66 




< 


1.19 


< 0.3 




Hydra-A 


0, 


.052 




3, 


,59 


± 


0, 


,10 


2, 


.15 


± 


0, 


,90 


0, 


,28 


±0.24 


3.7 


9.2 ±9.5 


A1060 


0. 


.011 


III 


3. 


,25 


± 


0, 


.05 














< 


0.04 


< 0.1 




A539 


0, 


.027 


III 


3. 


,21 


± 


0, 


.08 














< 


0.19 


< 0.1 




3A0335-h096 


0, 


.035 




3, 


,03 


± 


0, 


.07 


1, 


.42 


± 


0, 


,07 


0, 


,65 


±0.09 


8.0 


12.6 ± 1.9 


Virgo 


0. 


.004 


III 


2. 


,58 


± 


0, 


.03 


1, 


.25 


± 


0, 


,02 


1, 


,44 


±0.12 


0.5 


25.2 ± 2.3 


AWM4 


0. 


.031 


I 


2. 


,38 


± 


0, 


.17 


1, 


,14 


± 


0, 


,23 




< 


0.22 


< 0.2 




A400 


0. 


.023 


II-III 


2. 


,33 


± 


0, 


.14 














< 


0.13 


< 0.1 




A262 


0, 


.016 


III 


2. 


,20 


± 


0, 


.04 


1, 


.11 


± 


0, 


.07 


0, 


.25 


±0.19 


0.3 


6.3 ±4.2 



* The Bautz-Morgan morphological type. 

^ Hot-phase temperature and cool-phase temperature, both in keV. 

•f Ratio of mission integrals for the hot (sufBx h) and cool (sufRx c) phases, calculated over a central region of 
~ 100 h'^Q kpc, in unit of cm~'^. 

^ The 0.5-3 keV cool-phase luminosity in 10^^ ergs cm~^ s~^. 

* The volume filling factor of the cool phase in the central ~ 100 h^Q kpc, in units of percent. 
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Fig. 1 . . Ratio of the emission integrals between the cool and hot 
ICM components, calculated within ~ 100 h^Q of the cluster 
center and plotted as a function of the hot-component tem- 
perature. Filled and open symbols specify cD and non-cD 
clusters, respectively, as defined in subsection 2.3. Data refer 
to Table 1. 
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Fig. 2.. The cool component mass Mc of cD clusters, plotted 
against the 0.5-.3 kcV cool component luminosity Lc. Data re- 
fer to Ikcbc ct al. (1999) for the Centaurus cluster, Matsumoto 
ct al. (1996) for Virgo/M87, Xu et al. (1998) for AbcU 1795, 
Ikebe et al. (1997b) for the Hydra-A cluster, Ikebe (1995) for 
Fornax/NGC 1399, and Ikebe et al. (1997a) for Abell 496, 
Klemola 44, and MKW3s. The values of Lc employed here 
differ to some extent from those given in Table 1, mainly due 
to differences in the data integration radius. The data for 
non-cD ellipticals (filled circles) are from Matsushita (1997). 
The dashed line represents the scaling for a 1.0 keV plasma 
with a 0.5 solar abundance, assumed to have a constant den- 
sity of 2.5 X 10~^ cm""* with variable occupation volumes. 
Solid lines indicate the case when the plasma has a uniform 
but variable density and is confined within a constant spheri- 
cal volume of radius 25 kpc, with a filling factor of 1 (bottom), 
0.3 (middle), and 0.1 (top). 




Fig. 3.. Comparison between cooling flow rates measured with 
ASCA, against those reported previously (Edge et al. 1992). 
The ASCA values of M were taken from Ikebe et al. (1997b) 
for Hydra-A, Ikebe et al. (1999) for Centaurus, and Xu et 

al. (1998) for AbcU 1795; other ASCA values were calculated 
by substituting Lc in Table 1, after bolometric correction, 
into equation (8). Although ASCA errors arc shown only for 
Abell 1795 and the Hydra-A cluster, those for the remaining 
objects are similar. 
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Fig. 4.. The 0.3-5 keV cool-component luminosity Lc for a sam- 
ple of cD clusters measured with ASCA (taken from Table 1), 
shown against their central excess luminosity L^. The value 
of Lo refers to the integrated 0.5-3 keV luminosity of the nar- 
rower /3-component found by Mohr et al. (1999) through their 
double-/? fitting to the ROSAT PSPC profile; the integration 
is performed up to the maximum radius where the data exist, 
but the results are not very sensitive to the upper boundary. 
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Fig. 5.. Summary of the cases in which the hierarchical grav- 
itational potential has been observed; the Fornax cluster (for 
a distance of 20 Mpc; Ikebe et al. 1996), NGC 4636 (for a 
distance of 17 Mpc; Matsushita et al. 1998), the Centaurus 
cluster (Ikebe et al. 1999), Abell 1795 (Xu et al. 1998), and 
Abell 2199 (Xu 1998). For the latter three objects, = 75 
km s~^ Mpc~^ is assumed, (a) The integrated curves of the 
total gravitating mass. Abscissa, the 3-dimensional radius, is 
normalized to the interface radius Rip (see text subsection 3.3 
for definition), while ordinate is normalized to the mass con- 
tained within -Rip. Errors are not shown for clarity, (b) Re- 
lation between Rip and the total gravitating mass contained 
within it. These quantities are used to normalize the mass 
curves in panel (a). The straight line indicates a prediction 
for the universal halo (see text for detail). 
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scaled 3D radius 

Fig. 6.. Schematic illustration of three typical gravitational 
potential models for clusters, shown against the rescaled 
three-dimensional radius x. The King- type model has the core 
radius corresponding to a; = 1. The hierarchical double-/? po- 
tential assumes that the ICM emissivity is given by a sum of 
two /3-model components as expressed by equation (15), with 
the two core radii corresponding to a; = 1.0 and x = 0.12; the 
two components are assumed to have the same temperature 
and the same fj = 2/3, and the narrower emissivity compo- 
nent has a normalization 40 times as high as that of the wider 
emissivity component. This simulates the Fornax cluster. The 
Navarro-Frenk- White potential refers to equation (16), with 
f = 1 and Rs corresponding to x = 1.8; its normalization has 
been rescaled and offset to match the other two profiles over 
a range of x = 1 — 5. The models are no longer accurate for 
X > 5. 
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Fig. 8.. The ICM iron abundance in the cluster core region, 
plotted as a function of the cool-to-hot emission integral ratio 
there (Table 1). This updates Tamura et al. (1997). 
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Fig. 7.. The IMLR as a function of the system mass. Data for 
clusters and galaxy groups are taken from Fukazawa (1997), 
where the IMLR is calculated within a radius where the ICM 
density falls below 3 X 10~^ cm~^. Those for elliptical galaxies 
refer to Matsushita (1997), where the IMLR is calculated with 
4 times the optical effective radius. 
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Fig. 9.. Radial profiles of the IMLR (iron mass to light ration) for 
the Centaurus cluster (dashed curve, from Ikebe et al. 1999) 
and Abell 1060 (solid curve, from Tamura et al. 2000). The 
calculation properly takes into account the angular response 
of ASCA. 



